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An exciton theory for quasi-one dimensional organic materials is developed in the framework 
of the Su-Schrieffer-Heeger Hamiltonian augmented by short range extended Hubbard interactions. 
Within a strong electron-electron correlation approximation, the exciton properties are extensively 
studied. Using scattering theory, we analytically obtain the exciton energy and wavefunction and 
derive a criterion for the existence of a Bu exciton. We also systematically investigate the effect of 
impurities on the coherent motion of an exciton. The coherence is measured by a suitably defined 
electron-hole correlation function. It is shown that, for impurities with an on-site potential, a 
crossover behavior will occur if the impurity strength is comparable to the bandwidth of the exciton, 
corresponding to exciton localization. For a charged impurity with a spatially extended potential, 
CO ' in addition to localization the exciton will dissociate into an uncorrelated electron-hole pair when 

04 , the impurity is sufficiently strong to overcome the Coulomb interaction which binds the electron- 

hole pair. Interchain coupling effects are also discussed by considering two polymer chains coupled 
through nearest-neighbor interchain hopping t± and interchain Coulomb interaction V± . Within the 
t matrix scattering formalism, for every center-of-mass momentum, we find two poles determined 
^ ' only by V±, which correspond to the interchain excitons, and four poles only involving intrachain 

Coulomb V, which are intrachain excitons. The interchain exciton wavefunction is analyzed in terms 

■ of inter- and intra-chain character. Finally, the exciton state is used to study the charge transfer 
from a polymer chain to an adjacent dopant molecule. From a variational wave function for the 
total system, we explore the dependence of the probability of charge transfer on the acceptor level, 

' O ■ the hopping, and the wavefunction of the exciton. 

O ' PACS numbers: 71.35.+Z, 71.45. Cm, 71.55.-i, 71.27.-Ha 

> . I. INTRODUCTION 

■ In recent years, the exciton concept of electron-hole bound states has gained popularity in conjugated organic 
Qv^ polymers. Experimentally, it has been discovered that poly(phenylene vinylene) (PPV) and its derivatives, e.g., 

■ poly[2-methoxy, 5-(2' ethyl-hexoxy)-l,4 phenylene iiinylene] (MEH-PPV) can be used as the active luminescent layer 
f — I in electroluminescent light-emitting diode devices.EJ By using different conjugated polymers, pjolsjmer light-emitting 
0^ t diodes (PLED) have been fabricated which emit throughout the visible region of the spectrum.EJ"Q It is believed that 
■4-^ ' radiative recombination of excitons gives rise to luminescence, so a comprehensive understanding of exciton properties 
^ in polymer chains is very important to guide improvements in quantum efficiency of these PLED devices. Also, 
S conjugated poljaners have shown potential application in photonics for their large optical nonlinearity and ultrafast 

K [ response time.Q'B In view of the significant role excitons play in optical properties of a system, a detailed study of the 

(— ( ' exciton is also an important issue for polymer photonic device design. Theoretically, excitons in conjugated polymers 

Q ] are both attractive and challenging because of the coexistence of low-dimensional confinement and strong electron- 

O ■ electron (e-e) correlatiop. in. these systems. These ingredients have in fact led to many controversies during the study 

■ of excitons in polymers. ETlj j— ■ 
• ^ , Discussion of excitons in solids can be traced back more than half a century to the pioneering work by Frenkelt^l 

■ which even preceded than the band theory in solids. After numerous studies on excitons over several decades, the 
^ , exciton, as an elementary excitation, l^sJpeen well established in bulk insulators. This ^^he reason why in the 

■ - - ' existing exciton theories for polymers, lTEZI the standard exciton theory in semiconductorsEEl was usually borrowed 
with limited justification. In these theories, the polymer is regarded as a Peierls insulator which means the gap between 
the coruiuction and valence bands arises from the Peierls dimerization, a well-known nesting effect in one-dimensional 
metals-Ea Then the exciton states are solved considering the e-e interaction as a perturbation as in traditional exciton 
theory. But the polymer is significantly different from the conventional semiconductor. In a semiconductor, the 
electron correlation effects can largely be neglected and it is plausible to treat the e-e interaction as a perturbation. 
However, the polymer is typically a strongly correlated system with a moderate to large-polite Hubbard repulsion, 
and much of the band gap is due to the electron correlation rather than the dimerization.tJtj Thus the foundation of 
existing exciton theories in conjugated polymers is not so firm, and these theories have already led to some qualitatively 
incorrect results. For example, from these theories, the exciton is lower than the 2Ag, but typically the order 
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should be reversed. Again, the threshold of the conduction band is independent of U, the on-site Hubbard repulsion.tZI 
This is unreasonable since, when an electron is excited to the conduction band, double occupation must occur and 
cause an additional Hubbard repulsion. 

Having appreciated the central importance of electron correlatioitip polymers, several efforts have been undertaken 
to take account of the strong correlation effect on exciton states.cil E3 In these works, numerical exact diagonalization 
and density matrix renormalization group (DMRG) approaches are employed to handle small systems. From those 
numerical results, some useful information on the electron structure of the polymer can be captured by extrapolating 
the results to a long chain. However, these works cannot be regarded as constituting an exciton theory, since only 
some specific states are focused on in the calculations and the exciton is not treated as a quasi-particle. We must 
seek a more complete exciton theory of conjugated polymers, in which the correlation is stressed and the exciton is an 
elementary excitation. This kind of theory will then be useful to study the optical and transport properties related 
to the exciton. In this paper, we will develop a new exciton theory in the limit in which the Hubbard U is the main 
origin of the band gap (i.e., the polymer is regarded as a Mott insulator). In this regime, the spin has little effect on 
the exciton states and the energy difference between the singlet and triplet states is negligibly small compared with 
their binding energies. Using scattering theory, we will analytically calculate the exciton states and find a critical 
strength of the e-e interaction for the existence of bound exciton states. 

Currently, chemically synthesized polymers cannot be free from impurities and the "pristine" samples of polymers 
contain a non-negligible density of impurities and defects from the cross-linking, complex morphological effects, 
conjugation length effects, and some extrinsic defects. Such impurities sometimes critically influence the properties of 
the system, e.g., the transport. For example, the conductivity in traws-polyacetylene can be dramatically enhanced 
by 13 orders of magnitude by doping. Q In an impurity-free system, the exciton states with different center-of-mass 
momenta form exciton bands, and within the band the exciton moves coherently as a composite particle. The disorder 
tends to produce localized-states, and in one-dimensional systems, any nonvanishing impurity potential will lead to a 
localized electronic staterJ Thus it is interesting to examine the localization of an exciton, a composite particle, due 
to the impurities, and furthermore to determine if the exciton ceases to be a composite particle of electron and hole 
when the impurity is strong enough. In this paper, we will address the interplay of coherent motion of the exciton 
and different types of impurities in our conjugated polymer model. 

Strictly speaking, the polymer is only a quasi-one dimensional system, in which interchain couplings always exist, 
and sometimes their effects are striking. Many calculations have shown thaLjannlinear excitations like solitons and 
polarons may be unstable by taking the interchain coupling into account .EIEj Recently, great attention has been 
paid to the interchain effects in luminescent polymers, since manyj-experiments demonstrated that a large fraction of 
primary photoexcitations are interchain excitons or polaron pairsH3~Lj Also some theoretical calculations have been 
carried out to explore the interchain coupling effects on exciton properties .BEj Although the terminology of interchain 
exciton and intrachain exciton are widely used in current literature, these concepts are not so clearly delineated. From 
the principles of quantum mechanics, the wavefunction of every eigenstate in the coupled system must be distributed 
over the whole system, so it is difficult to distinguish from the wavefunction which one corresponds to the interchain 
and which one to the intrachain exciton. We will clarify what the interchain and intrachain exciton states are, and 
calculate their energies and wavefunctions. 

Photoconductivity in conjugated polymers can sometimes be greatly enhanced by intercalating or doping the poly- 
mer with a particular species.|Qf molecule. Interesting examples of this phenomenon occur when MEH-PPV is doped 
by fuUerene Ceo molecules. lZIlEI This is because the exciton, the bound electron and hole state in the polymer, will 
decay when the dopant molecule is introduced. The electron (or hole) in the exciton will transfer from the polymer 
chain to the doped molecule, giving rise to Bjicee carrier. Rice and Gartstein recently proposed a theory to explain the 
ultrafast time scale for this charge transfer .tj From a quantum mechanics perspective, assuming we have an exciton 
state in the polymer chain due to photoexcitation, when the coupling between the chain and the molecule is switched 
on, the electron will move in the whole system (including the chain and the molecule), and this state must have a 
lower energy than the initial state. So another point of view from which to study the charge transfer is to ask what 
percentage of the electron (hole) has transferred from the polymer chain to the dopant. This percentage should depend 
on the acceptor level and the coupling between the chain and the molecule, as well as the initial exciton wavefunction 
in the polymer chain. We will discuss this issue here. 

The paper is organized as follows. First, we develop an exciton theory for conjugated polymers in the strong 
correlation (large Hubbard U) approximation in Sec. II. In Sec. II. A we simplify a Peierls-extended Hubbard model 
to a model represented by spinless fermions with short-range e-e interactions in real space. Then we use t matrix 
scattering theory to determine the wavefunction and binding energy of exciton states analytically and derive a criterion 
for the existence of the exciton in Sec. II. B. This criterion is further proved according to the Levinson's theorem 
in scattering theory in Sec. II. C. In Sec. II. D a more formal and compact formalism for optical absorption in 
conjugated polymers is presented based on our exciton theory in the large-C/ limit. Sec. HI is devoted to the impurity 
effects on the coherent motion of the exciton. Using a suitably defined electron-hole correlation function, we study 
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different types of impurity. In Sec. IV, we investigate intercliain coupling effects by considering a two-chain system 
supplemented by nearest-neighbor interchain hopping t± and interchain e-e interaction V±. Using t matrix formalism, 
we analytically determine the poles corresponding to intrachain and interchain excitons, respectively. We also show 
the wavefunction of the interchain exciton. In Sec. V, the static Ag and Bu excitons are used to study the charge 
transfer in a molccularly-doped polymer. By constructing a variational wavefunction for the whole system, the energy 
of this variational state, and accordingly the probability of charge transfer, can be obtained. Finally, we summarize 
our results in Sec. VI. 



II. A NEW EXCITON THEORY IN CONJUGATED POLYMERS 

In existing theories, the polymer is regarded as a Peierls insulator, and then the exciton state is determined by 
treating the e-e interaction (including the on-site Hubbard interaction) as a perturbation. In this picture, the single 
particles (electron and hole) are defined based on a non-interacting Su-Schrieffer-Heeger model, so the band gap, from 
these theories, is independent of the e-e interaction. In a strongly correlated system, the electronic states are quite 
different from the non-interacting model. Since the ground state is half-filled, an electron excited to the conduction 
band must cost the additional Hubbard repulsion energy caused by the double occupation. In conjugated polymers 
and related organic conductors, it is now accepted that the origin of the band gap comes typically from the Hubbard 
repulsion rather than the Peierls dimerization. So the Hubbard term should be given priority when one develops an 
exciton theory. In this section, we will regard the polymer as a Mott insulator and develop the exciton theory in 
large-/7 limit. In this energy regime, double and higher order electron-hole excitations can be neglected because of 
their high energies (> 2[/), and a single configuration interaction approximation is reasonable in determining exciton 
states. Before carrying out the calculation, let us recall how large the Hubbard U (in units of electronic hopping 
energji,t)-|is in real systems: U/t ^ 3 — 4 in conjugated polymers and J7/t ~ 8 — 10 in segregated stack charge transfer 
salts.E3'cll Strictly, this strong correlation limit is applicable only when U ^ t, thus real conjugated polymers only 
marginally satisfy this approximation. 

A. Hamiltonians 

The theoretical model we consider is the Peierls-extended Hubbard model, i.e., the Su-SchriefFer-Heegei0 model 
augmented by an extended Hubbard interaction. For a one-dimensional chain, this model Hamiltonian is 

H = - {-iy5]{clci+,, + H.c.) + uY,nnni^ +vY,iPi - - 1) • (2-1) 

la I I 

Here cj^ creates an electron of spin a on site I, t is the one-electron hopping integral, ^ is a bond-alternation parameter, 

U and V are respectively the on-site and nearest-neighbor Coulomb interaction, nia- — cj^CAcr is the number operator, 
and pi = niij + n^. Since we will be concerned only with electronic excitations in this work, we consider a rigid 
dimerized ground state as a starting point, and do not specify its explicit origin (e.g., e-e interactions, electron- 
phonon couplings, or crystal structure). Strictly speaking, this Peierls-extended Hubbard model is directly applicable 
only to trans-polyacetylene. However, recent calculations have shown that the primary excitation in luminescent 
polymers like PPV can also be described within an effective linear chain model.E3 In these luminescent polymers, 
the lowest excitonic wave function extends over several repeat units. The properties of exciton are therefore not 
very sensitive to the detailed structure within the unit cell. From the viewpoint of renormalization, we can map 
the complex structure of a luminescent polymer into an effective Peierls-extended Hubbard system with the same 
significant physical properties by integrating out the superfluous degrees of freedom caused by the complicated unit 
cell structure. We have also neglected lattice relaxation, since many experiments and jtheories have demonstrated that 
e-e interactions dominate electron-lattice interactions in many luminescent polymers .e3 l3 This simplification enables 
us to handle e-e interactions in long chains and arrive at an understanding of electronic states in conjugated polymers 
without loss of essential physics, although the quantitative explanation of srnne lattice property like vibronic structure 
or bond length should, indeed, take into account lattice relaxation effects.E3 
To emphasize the electron correlation, we begin with the Hubbard model 

^o = -^E(cl'^'+i-+H-^-) + f^E"'T"U • (2-2) 

la I 

Altlmngh the exact wave function and the ground-state energy of this Hamiltonian have been obtained by Lieb and 
Wu,cZl the Green's function and correlation functions are difficult to calculate directly by using the exact wave function. 
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and it is also difficult to use their solution to study the exciton. As a practical alternative, here we make the strong 
correlation (large-C/) approximation. In this approximation, as we will see later--4he band gap is essentially C/, which 
is not the same as the exact solution for the Hubbard model by Lieb and Wu^ where the charge excitation gap is 
U — 4:t + ^'^^i{—^)"'[^nU — (t^ + |n^J7^)^/^]. However, for the strong correlation limit U '3> ^t, this difference is not 
important, and does not affect the exciton trends we wish to establish. 

The density product ni^nii can be expressed by the on-site electron number and spin operator, and the Hubbard 
model is rewritten as 

ifo = -t^(cLci+i. + H.c.) + [/^[ip,-(-iyS,.n,] , (2.3) 

la I 

where S/ is the electron spin operator at site /, 

S/ = ]^'^c\^(r„„,ci„, , (2.4) 

aa' 

and n; is a unit vector along the spin polarization axis of the electronS 

When U ^ t, the ground state is expected to have Neel order. We make an approximation by assuming that n; 
always coincides with the z axis, which implies that the spin excitations have been ignored. Thus, in this approxima- 
tion, the singlet and triplet excitons have the same energy and are not distinguishable. This is reasonable since for a 
one-dimensional Hubbard model, the spin and charge excitations are separated when [/ — s- oo and by order U ji for 
U ^ t, and the exciton is a charge excitation. This point is also directly justified-by comparing energies of the singlet 
and triplet states obtained from a finite chain exact diagonalization calculation.E2l 

Ho -tY,{4a^l+la+il-C.) + f - (-l)''^]cLc/^ ■ (2-5) 

la la 

This Hamiltonian is readily diagonalized by introducing 

Cla = -^y^ e*"[(Mfc -I- (TpiVk)aka + Pi{uk ~ (^PlVk)Pk] ■ (2.6) 
k 

Here the prime means that the summation runs over the reduced Brillouin zone |fc| < 7r/2 and Pi = (-1)'. Then 

ffo = -E' [(^^^ - ^)"L«fe- + (Ek + ^)Pljka] , (2.7) 

ka 

with 
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Ek = \l^+el, (2.8) 
£fe = -2tcosA:. (2.9) 



The functions Uk and Vk are 



Uk 



'I'k = ^1/1 




Ek V2 



In the case of J7 3> i,|-a| localized picture is more convenient. Two spinless fermions can be defined in the lattice 
representation as foUowscS 

ai=Y. 0{pi<j)VVNY,'e''''aka , (2.12) 

a k 

(3i ^Y.^{pia)^/2/NY^e''^'Pka , (2.13) 
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where 9{x) is the step function. Expanded in powers of t/U, cia- can be expressed by ai and /?/: 



t 



t 



cia = 9{pia)ai + 6{-pia)piPi + 6{-~pia) — {ai+i + ai-i) + 9{pia)pi — {f3i+i + fii-i) + 0(772) ■ 



If we include the bond ahernation part in our unperturbed Hamihonian, then 
H[,=Ho + Y,i-^ySt{clci+i^ + H.c.) 

la 

= jY.{h\hi + pIPi) + uY^plPi + ^ Y.^h\^2hi + P\Pi+2 + H.c.) 
/ I I 

+ StJ2i-hiPi+i + hi+iPi - Pj_^^h] + Pjh]^,) . 
I 

Here we have introduced the hole operator /ij = ai, and J — 2t^ /U . 
By introducing the Fourier transformations 



we can rewrite H'^ in momentum space as 

H[) = Yy^-^ ^ Jcos2k)hlhk + [U + J + J cos2k)pll3k + 2i5tsmk{hkP-k - /^Ifc^I)] 

k 

Making the Bogoliubov transformation 

r]k = cosdkhk +isin9kP^_k , 
7lfc = -isinOkhk + cos 6*^/31 j. , 



Hamiltonian (2.18) can be diagonalized if the relation 



sm29k = 



—4dt sinfc 
U 



is satisfied, yielding 



^0 = J2^'^''^kVk + tklllk) 



(2.14) 



(2.15) 



(2.16) 
(2.17) 

(2.18) 



(2.19) 
(2.20) 



(2.21) 



(2.22) 



with 



Ck = J(l + 5^) + ^(1 - 5^) cos 2k , 

lk = U + J{1- 3^2) + J(l + 3(5^) cos 2k . 



(2.23) 
(2.24) 



Operators r/^ and jj, create the hole and electron in the new valence and conduction band, respectively. Their lattice 
representations are 



li 



1 ikl 



rik 



(2.25) 
(2.26) 



which can be expressed by hi and j3i to order \/U: 
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St 

ll ^aj + ^jjihi+i-h^-i) . (2.28) 



r^i^hi~"-{pl,-pl,) , (2.27) 
St. 



u 

The inter-site Coulomb interaction is necessary to bind the electron and hole. We consider the ^-term in the 
Peierls-extended Hubbard model as a scattering potential, which has the local representation 

= vY,ih]+ih]hihi+, + /3;Vi/3/A/3z+i - hl^,pjPihi+, - Pl^^h\hiPi+{) . (2.29) 

Since the main interest here is an exciton, only the interaction between the electron and hole is relevant. To order 
we have 

Ht-t"" = -vJ2{vl,ilim+i + i\+Amii+i) ■ (2.30) 



B. Exciton states: t matrix theory 



Since the Hamiltonian is invariant with respect to translation, the exciton states can be classified according to the 
total quasimomentum K . We can write the exciton wave function as 



(2.31) 



where K is the center-of-mass momentum. The basis is chosen as 



(2.32) 



representing a created electron-hole pair from the ground state \g) with a separation *-in real space. We will determine 
the exciton state by using t matrix scattering theory. According to t matrix theoryt2l 



T{z)=u+ug{z)r{z) 



(2.33) 



where T{z) is the t matrix, Q{z) the notation for resolvent l/{z~ Hq), and U the potential operator. Equation ( 2.33| ) 
has the formal solution 



T{z)=U/[l-g{z)ll] 



(2.34) 



Using the basis of Eq. ( 2.32 ), we obtain the Green's function 
{r\g{z)\s) = G{r~s;z) 



s,K} — 



N ^ 



N ^ Z - {£k + e-k+K) 



(2.35) 



Here z — Ek + iO+ and the potential matrix is 



{s\U\s')^{^sMH':nt\^s'.K) 

= -VSss'{Ss,-l+5ss) 



(2.36) 



The utility of Eq. (2.34) rests on the possibility of actually constructing the inverse operator 1/(1 — QU). This can 
be achieved exactly in our case since, conveniently, the potential is of short range in the local representation. Actually, 
the portion of the potential U containing nonzero elements forms a 2 x 2 submatrix vmder the basis Eq. ( 2.32| ), 
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= -V 



1 
1 



Denoting D as the determinant of 1 — QU, we have 

D{Ek) 



1 + G{Q-Ek)V G{-2-Ek)V 
G{2-Ek)V 1 + G{0-Ek)V 



(2.37) 



(2.38) 



The determinant will vanish for some specific values of the energy, which are the energies of the localized states. 
Consequently, to find the energy Ek of the bound exciton state we look for the root of 



D{Ek) = . 

Subsequently, the wave function is calculated by solving the equation 

BrM = Y.^r\G{EK)\s){s\U\t)Bt,K 



(2.39) 



(2.40) 



First, let us focus on the static exciton, i.e., K — Q. In this case, the system is symmetric with respect to spatial 
inversion. Introducing the transformation 



BT = '^(^i ^ ^-') ' 

where Bi = Bi^k=o, and noticing 

G{s-t-Eo)^G{t-s-Eo) , 
we can write the determinant D as the product of two parts: 

D[Ea) = D^[Ea)D+{Eo) . 

Here 

D„ (i?o) = 1 + [G(0; Eo) - G(2; Eo)]V 
is for the Ag state with the wave function 



(2.41) 
(2.42) 

(2.43) 
(2.44) 
(2.45) 
(2.46) 
(2.47) 
(2.48) 
(2.49) 
(2.50) 

When 5 — Q, Eq does not equal the exact result by Lieb and Wr@. However, this is not a problem for estimating 
correct binding energies of the exciton states, since we will directly calculate the binding energy rather than first 
calculating the energy of the exciton. The Green's functions are readily calculated; for x > 



B- = ^V[G{1 ~ 1; So) - G{1 + 1; Eo)]B^ 



and 



D+{E^) = 1 + [G(0; So) + G(2; Eo)]V 
is for the B^ state with the wave function 

B+ = -V[G{1 - 1; So) + G{1 + 1; Eo)]B+ 
We denote x as the exciton binding energy: 

Eg - Eq , 

where Eq is the band gap: 

SG = (efc + efc)U=V2 = t^-4<5V. 



G(2Z + l;So) = , 



(2.51) 
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with 



G{2l;Eo)^^ ^ / i VI ^ \ ^2.52) 



- < 1 . (2.53) 



Then, by requiring D_{Eq) = 0, we obtain the binding energy of the Ag exciton as 

y2 



V + J{l + 5^) ' 



with corresponding wave function 



V-J{l + S^) 



and the wave function is 



(2.54) 



B- G{2l;Eo)-G{2l-2:Eo) 

v/-2[G'(0;i?o)"G'(2;i?o)] ' 

(2.56) 

Here G'{n;E) = ^G(n-E). 

For the i?„ state, when V > 2 J(l + (5^), the binding energy of the exciton is 

[T/-2J(l + 52)]2 



(2.57) 



^+ ^ G(2?;i;o) + G(2/-2;go) 

V-2[G'(0;i?o) + G'(2;i?o)] ' 
i?+ = . (2.59) 

When y < 2J(1 + (5^), a solution Eq < Eq satisfying equation D+{Eq) = cannot be found, i.e., there is no 
bound state. This result gives a criterion for the existence of a exciton. Finite chain numerical DMRG and exact 
diagonalization calculations show that the binding energy of the tS„ state is near zero when V is not large, but 
are not conclusive as to whether the state is strictly bound or free.Ell Our analytical results clearly indicate that the 
stable IBu exciton does not exist when V is less than 2J(1 + 5^). Since this critical value is half of the width of 
the continuum band, 4J(1 + (5^), and the bandwidth describes the kinetic energy of a free particle, this criterion is a 
reflection of the competition between the kinetic energy and the attraction of the electron and hole. 



C. Proof of the criterion: Levinson's theorem 

The criterion derived above can be proven by Levinson's theorem.El Namely, the number of bound states in repre- 
sentation s lying either above or below the continuum band can be determined by using 

5s{E,)~5,{Ef)=Tm, . (2.60) 

Here Ei is the lowest energy in the band and Ef the highest. Ug is number of states in any row of representation s 
separated £mm the band. 5a is the phase shift that appears in the usual partial wave expansion for the scattering 
amplitudejES which can also be extracted from the subdeterminants of det(l — QU), 

UnSa = -^Dj^Ds . (2.61) 

In the exciton case, 

E,=Eg , (2.62) 

Ef = (efe + h)k=o = Eg + 4J(1 + S^) . (2.63) 

It should be noted that in a one -dimensional system the density of states at the edge of the continuum band will 
diverge, so the form of Eq. (2.6C) in our case is 
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5,{E, - 0+) - Ss{Ef + 0+) = TTUs . (2.64) 



The Green's functions in different energy regions can be calculated according to the definition Eq. ( 2.35| ). We 



give the explicit expressions in Appendix A. Then the subdeterminants Z3_ and are obtained by straightforward 
calculations. 

For Eq = Ei — 0+, i.e., the energy is just below the onset of the continuum band, 

D_=-oo- iQ+ , (2.65) 



and the phase shifts are 



S-{E,-0+) =77 , (2.67) 

5+{E,-0+) = ^^ y>2J(l + <52). (2-68) 
For Eq — Ef + 0+, i.e., the energy is just above the top of the continuum band, 

D+ = +oo- iO+ , (2.70) 

and the corresponding phase shifts are 

S^{Ef +Q+)^0, (2.71) 

6+{Ef + 0+) = . (2.72) 

Thus we have 

6^iE,~0+)-S^{Ef+0+)^TT, (2.73) 
showing that there is always an Ag bound exciton for a nonvanishing V. However, when V < 2J(1 + 5^), 

S+iE,^0+)^S+iEf+0+) = 0, (2.74) 
indicating that there is no B.^^ bound exciton. The bound exciton will appear only when V > 2J(1 + 5^), since then 

6+{E^-0+)-d+{Ef + 0+)=TT. (2.75) 

To get an overall picture of the phase shift in this one-dimensional system, we also examine the phase shifts at 
Ei + 0+ and Ef — 0^ . For Eq — Ei + 0+, i.e., the energy is just above the onset of the continuum band, 

= l + (2.76) 



and the phase shifts are 



6-{E^+0+)^tt/2 , (2.78) 
S+{E. + ) = |^ V>2Jil + S'). (2.79) 



For Eq = Ef — 0+, i.e., the energy is just below the top of the continuum band, 
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D- = l+ ^ - «0+ , (2.80) 



2J(l + 52) 
V 

2J(l + (52) 



^+ = 1 - 7777rTT2T - ' (2-81) 



and the phase shifts are 



5-{Ef . (2.82) 

5+{Ef-Q ) -i^^l^ F>2J(l + 52). (2.83) 

Now let us study the behavior of the phase shift as a function of energy. For the Ag state, when the energy passes 
through the edge of the band the phase shift falls discontinuously from tt to 7r/2, then gradually approaches zero at 
the top of the band. So there is always a bound state. For the exciton, when V < 2J(1 + 6'^) the phase shift 
increases from zero as the energy increases from the bottom of the band, and approaches tt/2 just below the top of 
the band, and then drops to zero again when we pass through the top edge of the band. Thus no bound state exists. 
When V > 2J(1 + S"^), the phase shift decreases from tt to tt/2 as the energy increases from the bottom to the top of 
band, and abruptly falls to zero when we cross the edge of the band. Thus a bound state appears. The discontinuity 
of the phase shift at the band edges is due to the infinite density of states at the bottom and top of the band. 

Figure 1 shows the wave functions of the Ag and S„ states with U — IQt, V — 0.5t, and S — 0.2, corresponding 
to binding energies x+ = 0.024i and X- — 0.353i, respectively. We can see that the wave function of the Ag exciton 
decays more rapidly than that of the one. For x > 0, we introduce the parameter z, 



z = - ln(^ — j > . (2.84) 

Since for large I, 

G{2l;Eo) - e-^|2'l , (2.85) 



we can define the width R of the K = Q exciton by i? = 2/ z. From Eq. ( p.84|) , we estimate the width is about 3 
lattice constants for the Ag exciton and about 12 lattice constants for the as shown in Fig. 1. 

We calculate the energy of the exciton for if 7^ from D{Ek) =0. A straightforward computation of Eq. ( 2.35 ) 
gives 

G{0-Ek)^- J— ^ , (2.86) 
G{2-Ek) = G*{-2-Ek) , 

where 

a = a; + 3J + 5JV + J(l-(52)cos2iC , (2.88) 
h= J{1- 5'^)sm2K , (2.89) 
c = a; + J(l-(5^)- J(l-(5^)cos27f . (2.90) 



By solving the equations 



1 + G{Q-Ek)V +\G{2-Ek)\V , (2.91) 
1 + G{Q-Ek)V -\G{2-Ek)\V , (2.92) 



we can compute the energy of the moving excitons. In Fig. 2, we have described the energy of the excitons as a 
function of the center-of-mass momentum K . There are two branches in the energy spectra. The energy difference 
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between these two branches becomes smaller when K increases, and reaches a minimum at K = 7r/2. The bandwidths 
of these two branches are approximately J. 

When V ^ J, pur results show that both Ag and i?„ excitons have binding energy V, which is consistent with 
physical intuition.Ej In the strong coupling limit U ^ V '3> t, Guo et al. used a local (zero hopping limit) picture 
and argued, since the ground state has all sites singly occupied and the exciton states are linear combinations of 
configuration ...11120111..., where the numbers denote site occupancies, that the exciton energy \s U — V . The 
electron-hole continuum consists of all states in which the double occupancy (electron) and the empty site (hole) are 
separated by more than one site (e.g., ...11211. .1011...), which has energy U . Thus the binding energy is V . Another 



prediction from our theory is that the 2Ag st 
many non-luminescent conjugated polymers.! 
reason for this ordering in our model. 



IS a lower energy than the li?„, which is the observed ordering in 
3 The strong Coulomb interaction regime we consider here is the 



D. Optical absorption 

We will calculate the optical absorption from Fermi's golden rule 

a{Lo) cx - V|(n|J|5)|25(c^-^„) , (2.93) 

LO ^ — ^ 

n 

where J is the current opertor and \n) is the excited state with one electron- hole pair. This expression can be written 
in a more general form if we denote — J|(?), 

a(w) cx — 3_ lim 1 . (2.94) 

TTtj e^o+ Lu + ie — H 

Hamiltonian H referred to here is that determining the energy of the electron-hole pair with K — 0. The current 
operator in the polymers reads 

J - zi^[l - i-iyS]{cl,^ci, - 4q+i,) . (2.95) 
Using the spinless fermions rji and 7/ defined in Sec. II. A, we rewrite Eq. ( ^.95 ) as 

Thus the optical absorption can be expressed by the electron-hole Green's function 

a{uj) cx 3[G(0; lj) + G(2; lu)] , (2.97) 

where G{l\uo) is the Green's function of H, satisfying 

g{z) = g{z) + g{z)T{z)g{z) . (2.98) 

If we denote G{n;uj) — G„, then 

a{u:) cx -—^{Ga + G_2) - 7^[(Go + G2)' + V{GI - G^GI + G^Ga - G^)]} (2.99) 

and D{u:) = [1 + (Go - G2)t^][l -f (Go ^ Ga)^]. 

From Fig. 3, the exciton has acquired 52% oscillator strength when U = 5t and V = t. If we increase V and 
thus have an exciton with a larger binding energy, the i?„ exciton will gain more oscillator strength. For U ^ 5t 
and V = 2t, the i?„ exciton has achieved 95% strength, and the strength of the transition to the continuum is 
correspondingly diminished, as shown in Fig. 4. The large transition strength for the exciton state is a characteristic 
feature of one-dimension. 
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III. IMPURITIES AND THE COHERENT MOTION OF EXCITON 



As stated in the introduction, "pristine" samples of a polymer cannot eliminate all impurities and defects. Moreover, 
the fluctuatioas (both quantum and thermal) of the lattice are also a kind of intrinsic disorder for the electronic states 
in polymers.E^ Excitons represent a coherent composite particle motion of correlated electrons and holes, whereas 
impurities tend to produce localized wavefunction. We know that an impurity has a strong effect on transport 
properties in many systems, especially low-dimensional materials, so a natural question arises: how do the impurities 
affect the coherent motion of the exciton? l— . 

In polymers, two kinds of impurity are often referred to in the literature.Ej A site impurity is represented by a local 
potential at site 0: 

H,=VoY.clc,, , (3.1) 

CT 

and a bond impurity which acts on the bond between sites and 1: 

H2 = -Wo + H.c.) . (3.2) 

cr 

Both of these impurities have very localized (on-site) potentials. However, for a charged impurity, its potential, in 
principle, may be of long-range. Thus we have two length scales here: one is the range of the impurity potential (li) 
and the other the range of the (screened) Coulomb interaction (ly)- The latter is equal to one lattice constant in our 
model. Since the impurity competes with the Coulomb interaction (in exciton states) differently in the two regimes 
{Iv > li on Iv <li), impurity effects are expected to be different. 



A. On-site impurity potentials 



For the site and bond impurities described above, we can rewrite them in the spinless fermion representation, giving 

Hi^Vo{~vlvo+lho) . (3.3) 



H2 = -W^o(-»7o7i - llvl + llvl + mio) ■ (3.4) 

Since H2 involves the creation and annihilation of an electron-hole pair, it must be less important than Hi by order 
1/U. This can be seen more clearly by using the unitary transformation = e~^He^ , 

= £/ + 2J(l-^2) ("'^o^o - vim - 7^70 - ihi) • (3.5) 

Although the site impurity seems more realistic from the above analysis, we will study three kinds of impurity to 
arrive at a unified picture of impurity effects: 

hT" = Vo{-vlvo+jho), (3.6) 
Hr = Voi-vlvo-lho). (3.7) 
Ht" = yoivlrio+lho)- (3.8) 

Hamiltonian in which the impurity attracts the hole and repels the electron, or vice versa, imitates a local 

charged impurity. Hamiltonian -ff™'', in which the impurity attracts both the electron and hole, acts as a trap for 
particles. In Hamiltonian i?™'', the impurity potentials are repulsive for both the electron and hole, describing a 
barrier effect. The last two types of impurity can be viewed as simulating the cross-linking and conjugation breaking 
effects in conjugated polymers. 

There is no translation invariance once the impurity is included, so we will work in real space and the Hamiltonian 
we must study reads 
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+ <5^)(7/7i+2 + 7/+27O} +"^^^{4+1^1^1^1+1 + 7/+i7/7i7i+i 

I 

- vl+iihm+i - ij+ivjmii+i) + hI"""^ . (3.9) 

The key issue here is how to measure the coherence in the excitonic composite particle. We can do this by defining 
the correlation function between the electron and hole in the lowest state in the one electron and one hole subspace: 

where the deviations are 

SA = A-{A), (3.11) 
and the density operators of electron and hole are 

Phil) = vim , (3.12) 
Peil)^ljli- (3.13) 



In the impurity-free system, this correlation function ( jj.lO ) is connected with the relative wavefunction of the lowest 
exciton state |*o) = ^ Bsjj+.vlld) by 

7^f■^<=°(^,0 = l^i-H' • (3.14) 

When we add an impurity, the correlation is expected to decrease. The closer TZ is to 71^'^°'^, the more correlated 
are the electron and hole in the lowest excitonic states, while TZ approaching zero means that there is no correlation 
between the electron and hole; in other words, this excitonic state has lost all its coherence. 

The effects of the first kind of impurity are illustrated in Fig. 5, which shows the electron-hole correlation functions 
for different sites in a finite system of size N = 10. We emphasize that the parameters we use [U = IQt and V = t) 
ensure that the exciton has a very localized relative wavefunction, so that finite system corrections and the boundary 
condition effects are not important. In Fig. 5, all the correlation functions exhibit a crossover behavior around 
Vo ~ Q.lt. This can be understood in term of a relevant energy scale of the exciton, namely, the width of the exciton 
band, which equals J — 2t^ /U . This crossover behavior, which occurs at Vb ^ J, describes the localization of the 
exciton, i.e., the free exciton becomes trapped. We can also calculate the charge density at the impurity site as the 
impurity strength increases. Since in an impurity-free system the hole (electron) is uniformly distributed, and from 
our exciton theory the electron and hole do not tend to occupy the same site, the hole density at the impurity site 
is 2/N . We see a crossover again in Fig. 6 when Vb is comparable to the bandwidth J. After this crossover, the 
hole density at the impurity site approaches 1, clearly showing that the exciton is trapped by the impurity, and the 
correlation between the electron and hole gradually vanishes, although, as indicated in Fig. 7, they are bound together 
near the impurity. 

For the second type of impurity, from Fig. 8, a crossover is also observed if the impurity strength is similar to 
the exciton bandwidth, again indicating that the exciton is trapped. But when Vq is larger than V , the correlation 
function abruptly falls to zero, which implies the total breakdown of the exciton as a composite particle. This is 
because, when Vb is large enough, it is a lower energy for the impurity to trap the electron and hole separately rather 
than the impurity trapping the hole and then the electron being trapped near the hole due to the Coulomb interaction 
(as for the first kind of impurity). Thus the electron and the hole occupy the same site and they have no Coulomb 
interaction. This is not an exciton. 

Now we consider the third type of impurity. The correlation function behaviors in Fig. 9 seem more complicated 
than for the other two types of impurity. The correlation function in which the hole is at the impurity site shows 
an analogous crossover behavior when Vb is near the exciton bandwidth J to the first and second types of impurity. 
However, if both the electron and the hole are left (or right) of the impurity, from the correlation function we see that 
they have not felt the impurity. On the other hand, if the electron and the hole sit on different sides of the impurity, 
there is a crossover at Vb ~ J, but part of correlation between the electron and hole survives. 

The different effects of these three kinds of impurity can be further understood if we project the lowest excitonic 
state to the free exciton states with momentum K. We depict in Figs. 10, 11 and 12 the distribution \Zk\'^, where 
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Zk = {^k\^), and \'^) is the lowest excitonic state in these three disordered systems. In the impurity- free system, the 
lowest state is the linear combination of exciton states with K = and K = tt (they are degenerate) . In the presence 
of an impurity, the exciton state will be scattered to other exciton states with different K and the distribution of 
momenta will broaden from the S function. A free exciton state with a specific K can be defined only when the width 
of the distribution of momentum is narrow enough the lifetime of this state is long enough). This is analogous 
to the quasi-particle in Landau's Fermi Liquid theory.E3 From these figures, we see that for the first and second types 
of impurities, after the crossover at Vb ^ J the distribution in momentum space becomes so broad that we can hardly 
identify the original exciton state with momentum K = or K = n. For the second type of impurity, when Vq is larger 
than V, the final state has no distribution at all on any free exciton state, also indicating that the final state is no 
longer excitonic. However, for the third kind of impurity, after the crossover at Vq ~ J, the distribution in momentum 
space still has two sharp peaks at i^T = and K = n. This is the reason why the exciton is still coherent, as shown in 
the correlation functions, and in this sense the exciton can be regarded as a quasi-particle in this disordered system. 



B. Extended impurity potentials 

For a charged impurity, the range of (screened Coulomb) potential can be extended over several lattice constants. 
As an illustration, here we consider a specific impurity potential 

iJ^'°P = Vo{-iil7]o + 70 7o) + y(-'7I7/i + 7l7i - + 7li7-i) • (3-15) 

Its range (k) is three lattice constants which is longer than that of the Coulomb interaction {Iv)- The correlation 
functions are illustrated in Fig. 13, from which we observe again a crossover around Vq ^ J, indicating the free 
exciton becomes trapped. Interestingly, when Vo is sufficiently large compared to V, the correlation functions falls 
further abruptly to zero, indicating the dissociation of the exciton into an uncorrelated electron-hole pair. Note that 
this does not occur for the charged impurity H™^^ with the on-site potential, since the impurity range is then less 
than the trapped exciton size. From charge densities shown in Fig. 14, we find that for Vq = 0.5t (just after the first 
crossover), both the electron and the hole (thus the exciton) are trapped around the impurity. For Vq = 5t (after the 
correlation goes to zero), the hole is trapped by the impurity while the electron is repelled from the impurity. The 
dissociation of excitons here is easily understood. Because the impurity attracts the hole and repels the electron, when 
the impurity strength becomes sufficiently strong, the Coulomb attraction cannot overcome the impurity repulsion 
to bind the electron and hole together. This impurity-induced exciton dissociation may be invoked to interpret 
impurity-enhanced photoconductivity observed in certain experiments. We can project the lowest excitonic state in 
the system with the impurity to free exciton states with different momenta. In Fig. 15, the distribution in momentum 
space changes from a sharply localized one (Vq = O.lt, before the crossover) to a very broad Gaussian distribution 
(Vo = OM, after the crossover) and finally goes to zero (Vq = 5t). This is consistent with the picture that the free 
exciton becomes trapped, then dissociates into an uncorrelated electron-hole pair with increasing impurity strength. 



IV. INTERCHAIN COUPLING AND INTERCHAIN EXCITONS 



The interchain coupling can strongly infiuence the energy and stability of the nonlinear excitations such as solitons 
and polarons. Current interest in interchain coupling and the intrachain and interchain exciton crossover in polymers 
stems from the experimentally observed large amount of interchain excitations in luminescent polymers like PPV. 
However, the concept of an interchain exciton and how to distinguish the interchain and intrachain excitons, arc not 
very clear. The wavefunction is not so useful to specify whether a state is an interchain or intrachain exciton, because 
the wavefunction of any state, in principle, will spread over the whole system if interchain coupling is present. 

To demonstrate interchain exciton states in our approach, we study a two-chain system coupled by the nearest- 
neighbor hopping, 

i^hop = -t±J2'^c{i^^2ia + H.c.) , (4.1) 

la 

and nearest-neighbor interchain Coulomb interaction, 

Hcon = Vi^ Y.^pu ~ l)(p2i - 1) . (4.2) 
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Here 1, 2 are chain indices. Now the unperturbed Hamihonian in the large-[/ Umit becomes 
Ho ^ Hq + ifhop 

= J2^'^''^^^kVtk + hllklik) + ij- ^{vlm - l\ll2l + H.c.) . (4.3) 

In momentum space, this is 

^0 = ^{^kvlk^tk + hllkltk) + t± ^ivlkVik ~ llkl2k + H.c.) , (4.4) 

i,k k 

which can be readily diagonahzed, yielding 

Ho = J2(^kVlkVik + Eilllik) ■ (4.5) 

I,k 

Here / = 1, 2 is the band index (hereafter we use smaU i,j for the chain indices, and capital /, J for the band indices). 
The details of the diagonalization are given in Appendix B. Now we have two conduction and two valence bands with 
the dispersion relations 

El'^ = ek T t± , (4.6) 



eI'^ = ik T t± . (4.7) 



The interaction Hamiltonian now contains two parts 

Hv = H,nt + Hcon ■ (4.8) 

The electron- hole interaction part of iJy, which is relevant to the exciton state, is 

H'^^ = -V"^(77|;+i7|;7,i?7.i;+i + lll+llilr|^a^l+l) - Vi_ Y^{ri\illil2imi + llivlimaii) ■ (4-9) 



If we define local fermion operators 



'?^' = :^Ee''''^^^ ' (4.10) 

k 



we can rewrite Eq. ([4.9|) as 



2 

ij I 



^ Y^{v\i+ii\ii2imi+i + f}\i+iilm^^ +1-^2) 



2 



^ "fWjiVjnii - J2'^vuilii2imi + vliilnimi + 1-^2). (4.12) 



2 2 



Thus we can construct the exciton wavefunction in this two-chain system as 

Here the center-of-mass momentum K is still a good quantum number, and the basis is 

\i'l'K) = ^,Y.^''''l\i+s€i\9) ^ (4.14) 
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which means the exciton state is a combination of every possible electron-hole excitation in different conduction and 
valence bands. The free electron-hole pair Green's function under the basis Eq. ( 4.14 ) is 



1 



zk - H( 



= Sir5jj,G"{r- s;zk) 



(4.15) 



1 x-^ e''^' 



N^z^-{El + Ei,^^) 



and the scattering potential is written as 



V 

Y 



2 



(4.16) 



(4.17) 



with F{I, J; /, J) = 1 and F(l, 1; 2, 2) = F(2, 2; 1, 1) = F{1, 2; 2, 1) = F(2, 1; 1, 2) = 1, and otherwise F = 0. Wc can 
solve for the exciton states by locating the roots of the determinant det{l — QU) according to t matrix theory. The whole 
determinant can be decomposed into blocks by appropriate transformations and we achieve two subdeterminants, 



Di 



1 + -^G"(0;zk) ^G"(0;zk) 



V 



V, 



G''\0;zk) 1 + ^G^\0;zk) 



and 



Do 



1 + Y±g'^0-zk) ^G''{0;zk) 
y^G^\Q;zK) 1 + ^G^\Q;zk) 



(4.18) 



(4.19) 



Di and D2 are determined only by interchain Coulomb interaction Vj_ corresponding to interchain excitons, and two 
subdeterminants, 



and 



D4 = 



1 + ^ [G" (0; zk) + (0; zk)] ^ [G" (-2; z^) + G^' (-2; zk)] 
^ [G" (2; zk) + G22 (2; zk)] 1 + [G" (0; zk) + G^^ (0; zk)] 



1 + ^[G'\0;zk) + G^\0;zk)] -^[G12(-2;zk) + G21(-2;zk)] 
^ [G'^ (2; zk) + G^i (2; zk)] 1 + [G^^ (0; zk) + G^^ (0; zk)] 



(4.20) 



(4.21) 



which are determined only by the intrachain Coulomb interaction V corresponding to intrachain excitons. Equations 
Di = and D2 = have a single root, respectively, whereas both D3 = and D4 = have two roots. Thus there are 
a total of six exciton bands: two interchain exciton bands and four intrachain exciton bands. Figures 16 and 17 show 
these intrachain and interchain exciton bands. The relative energy ordering of the interchain and intrachain excitons 
depends on the ratio V±/V. 

It is interesting to study the wave function of the interchain excitons. The static interchain exciton can be repre- 
sented in real space as 



(4.22) 



tj s • ; 
For the lower exciton state determined by Di{zq) — 0, we obtain 

G"(s;^o)-G22(s;zo) 



All _ _J^l-Z _ _ 
-Tic. -^c 



2V-[Gii(0;^o) + G22(0;zo)]' 



(4.23) 
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which represents the ampHtude for the electron and hole being on the same chain; while 

2v/-[G"(0;2o) + G22(0;zo)]' 
is the amplitude that the electron and hole are on different chains. Here 

Gii(0; z) = ^ (4.25) 

6-22 (0; z) = ^ , (4.26) 

^{Eg - z + 2t^){EG -z + 2t^ + W) 



where Eg is defined as Eq. ( 2.50| ) and W = 4J(1 + S^). Figure 18 illustrates the intrachain wavefunction Al^ and 



interchain wavefunction A^^ for this interchain exciton state. We can see that although the state is for an interchain 
exciton, there is still some probability for the electron and hole to be on the same chain. 
For the higher interchain exciton determined by -D2(^o) = 0, we have 

All ^A22^Q^ (4 27) 

^21 ^ ^12 ^ G^^js; Zq) 

v/-[Gi2(0;zo)]' ' 

Gi2(0; z) = ^ . (4.29) 

^{Eg-z){Eg-z + W) 

In this interchain exciton, there is no relative amplitude between the electron and hole if they are on the same chain. 
In Fig. 19, we plot the wavefunction A"^^ for this interchain exciton. 

For more complicated Coulomb interactions within and between the chains, the exact interchain and intrachain 
exciton poles ar e diffi cult to obtain analytically. Instead, we can measure the correlation between the electron and 



hole using Eq. ( 3.10 ). In Fig. 20, we depict the intrachain and interchain electron-hole correlation functions for our 
simple interchain coupling situation for two iV = 12 chains. Here we choose a fixed interchain hopping t±/t = 0.2. 
The transition at V±/V =1.18 shows that the lowest exciton state changes from an intrachain exciton to an interchain 
one. 



V. CHARGE TRANSFER IN A MOLECULARLY-DOPED POLYMER 



Photoinduced charge transfer from a polymer chain to an adjacent dopant molecule, such as in PPV/Cgo blends 
has attracted much recent attention, because it can greajtiy increase the photoconductivity in polymers. In recent 
theoretical work on this phenomenon. Rice and GartsteinLJ proposed a mechanism to explain the observed ultrafast 
time scale of this process. In this section, instead of discussing the time scale, we attempt to calculate the final state 
wavefunction of the whole system comprising the polymer chain and dopant molecule. This can tell us what part of 
the electron in the exciton has transferred from the chain to the dopant. To make our idea more transparent, let us 
briefly describe this photoinduced charge transfer process. In the ground state, there is no overlap between the chain 
and the dopant. The photoexcitation produces an exciton state in the polymer chain. Then, due to the coupling 
between the polymer and the dopant molecule, the electron (or hole) will transfer from the chain to the adjacent 
molecule. As a simplified Hamiltonian, we consider 

H = Hchmn + Ae ^ cJ,jC„, + i/tran • (5.1) 
m 

Here we are modeling the dopant molecule by assuming it has an acceptor level with energy Ae , which couples to the 
polymer chain only by nearest-neighbor hopping 

Htr.n = ~vJ2{cIjo + H.c.) = ^(cL7fe + H.c.) . (5.2) 
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A schematic diagram of our model is shown in Fig. 21. The whole system consists of a polymer chain and Nd dilute 
noninteracting dopants. Before the coupling between the polymer chain and the dopant is switched on, the system 
has an exciton state on the polymer chain. When we turn on the coupling, the electron in the exciton will transfer 
between the chain and the molecule. Thus we can construct a variational wavefunction 

\^) = a|*o> + 7^ E E «fe47fc|*o) , (5.3) 

with the condition (^'|^') = 1. The first term describes the electron remaining on the chain as a component of the 
exciton and the second term describes that the electron with different momentum has different probability to transfer 
to the dopant molecule. Here j^I'o) is the assumed static exciton state within the polymer chain. 



(5.4) 



which can be represented by the relative momentum between the electron and hole 



Its energy is 



(«'o|i?chain|«'o) = Eo . (5.6) 

The variational state |^) must have a lower energy than Eq , 

e = (^liJl^-) - So < . (5.7) 

de de 

From — — — and — — — 0, we obtain two coupled equations: 
oa oak 

lv7 

-j^av ^ ak[e + Eo - {ek + A,)] , (5.8) 



= ^E«^l^'»l'yf- (5-9) 

k 



Using Eq. (5.8), can be eliminated from Eq. ( |5.9| ) and we have the eigenvalue equation for e 

e = c— > 7 ir^ — (Q ■ (5.10) 

N ^ e + Eo-i^k + Ae) ^ ' ^ ' 

Here c = N^/N is the dopant concentration. Once we have found the negative solution of e, the probability that the 
exciton remains on the chain is 



where F'{t) — ' . So the probability of charge transfer is 



P = l-a^= . (5.12) 

1 - cF'(e) ^ ' 

For the i?„ exciton in the polymer chain 

Bfe = E5i^e-*". (5.13) 

If we assume for demonstration purposes that the exciton is highly localized (this is not necessary in our theory), i.e., 
B+ = Bt^ = and = for other /, then 
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Bk = V2cosk, (5.14) 

and we can write F{e) in a very compact form: 



with 

E ^ -{e + Ea - Ae -2J5^) . (5.16) 

Figure 22 illustrates the probability of charge transfer to the dopant. The dopant concentration is set to be 0.2, 
a typical value for a molecularly-doped polymer. We see that when the acceptor level is near the exciton energy 
Eq, a crossover will occiu". When the acceptor level is below this crossover value, the electron is mainly on the 
dopant. Otherwise, the electron is mainly on the polymer chain. The coupling strength v affects the charge transfer 
by controlling the width of the crossover. The smaller v is, the more rapid is the crossover. 

For the Ag exciton, if we make the assumption = —BZi = then 

\Bk\ = \Y,Bfe''''\ = V2smk, (5.17) 
I 

and we have 



We plot the charge transfer probability for the Ag state in Fig. 23. We can see that there is a threshold for Ag. Below 
this value, the electron will thoroughly transfer to the dopant molecule. However for the state there is a long tail 
below the critical value, indicating some fraction of the electron can still be found in the polymer. Having gained the 
knowledge of how A^, v, and the exciton wavefunction influence the probability of charge transfer, one will be able 
to control this transfer process in the conjugated polymer. 



VI. CONCLUDING REMARKS 



In this paper, we have extensively studied the exciton states in conjugated polymers by emphasizing the dominant 
role of e-e correlations. The model we studied here is the widely- used Peierls-extcnded Hubbard model with frozen 
bond dimerization. First, in the lavge-U approximation, we mapped this model to a spinless fcrmion model with 
only nearest-neighbor Coulomb interaction in real space. The short range interaction enabled us to apply t matrix 
theory to analytically calculate the energy spectrum and wavefunction of bound (exciton) states. We have found 
that there always exists a stable Ag exciton as long as the nearest-neighbor Coulomb V is nonzero; for the Bu state, 
however, a stable exciton state can exist only when V is larger than the half width of the continuum band. This 
criterion has been proven based on Levinson's theorem. In our results, we have a correct ordering for Ag and i?„ 
states, i.e., 2Ag < IB^, as observed in most conjugated polymers. The impurity effects on the coherent motion of the 
excitons were also investigated in this large-C/ approximation. The coherence of the exciton can be measured by an 
appropriately defined electron-hole correlation function. We have studied impurities with on-site potentials as well as 
a charged impurity with a more extended potential. There are three kinds of impurity with the on-site potential: the 
first is like a local charge, attracting holes but repelling electrons (or vice versa); the second acts as a well, attracting 
both electron and hole; the third is like a barrier which repels both electron and hole. We have found that for the 
first and second type of impurities, the electron-hole correlations exhibit a crossover when the impurity strength Vq is 
comparable to the exciton bandwidth J, which describes the exciton being trapped by the impurity. For the second 
type of impurity, if the impurity strength is larger than the Coulomb interaction V, the deep well will trap the electron 
and hole separately, leading to the total de-correlation of the exciton as a particle. For the third type of impurity, the 
exciton coherence can survive the impurity and the distribution in momentum space has a sharp peak which means 
the exciton still moves freely. For the charged impurity with an extended potential of range greater than ly, the 
range of the Coulomb interaction, the free exciton becomes trapped at Vq ^ J, analogous to the situation for on-site 
impurity potentials. However, unlike the charged impurity with the on-site potential, the exciton dissociates into an 
uncorrected electron- hole pair when Vq is sufficiently large compared to the Coulomb strength V. 
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We have also investigated the effects of interchain couphng and the resulting interchain exciton states within the 
strong correlation approximation by considering a two-chain system with nearest-neighbor interchain hopping t± and 
Coulomb interaction V± . In this coupled system, we have two conduction bands and two valence bands. Within the t 
matrix formalism, we have found six poles for every center-of-mass momentum K , in which two poles are determined 
solely by V±, corresponding to interchain excitons, while the other four poles are determined solely by the intrachain 
Coulomb interaction V, corresponding to intrachain excitons. We have also illustrated the wavefunctions of the static 
interchain exciton. There is still some amplitude for the electron and hole being on the same chain for the interchain 
exciton state. For more complicated Coulomb potentials, we propose a way to distinguish the interchain and intrachain 
excitons, namely by comparing the interchain electron-hole correlation function with the intrachain one. 

The charge transfer in a molecularly-doped conjugated polymer has been studied by constructing a variational 
wavefunction for the whole system including the polymer chain and dopant molecule. We modeled this coupled 
system by simply regarding the molecule has an acceptor level which interacts with the polymer chain by nearest- 
neighbor hopping V. Minimizing the energy of the state, we have obtained the energy of the variational state and, 
accordingly, the probability of charge transfer. We have shown that a crossover behavior will occur when the acceptor 
level is near to the exciton energy. When the acceptor level is higher than this crossover value, the electron mainly 
remains on the polymer chain. Otherwise, most of the electron density will transfer to the dopant molecule. The 
hopping V controls the width of this crossover, the larger v is, the more gentle is the crossover. The wavefunction is 
also an important influence on the charge transfer. For the Ag state, there is a threshold for the acceptor level. If Ae 
is less than this value, the charge transfers to the molecule thoroughly and the percentage of electron density in the 
polymer chain is zero. 

Our calculations in this paper presented a comprehensive picture of the exciton in conjugated polymers, in a limit 
in which the electron correlation effects have been taken seriously and consistently. Our exciton theory can be readily 
extended to a system with a relatively long range Coulomb interaction. Also using our spinless fermion representation 
for the Peierls-extended Hubbard model, biexciton states can be obtained either by the Heitler-London method 
or diagonalization of the Hamiltonian in two electron-hole pair space. Although for real conjugated polymers, the 
Hubbard U is not so strong and our results cannot quantitatively match the energy levels in luminescent polymers, our 
theory is useful for understanding several puzzles which have arisen from correlation effects in conjugated polymers. 
Finally, we note that recent experimeatal evidence has demonstrated that there is an excitonic contribution to the 
pairing mechanism in YBa2Cu307_5E3 We expect that our exciton theory can give some guidance for exciton effects 
in high-Tc superconductors by extending the formalism to two dimensions. 
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APPENDIX A: EXPLICIT EXPRESSIONS OF GREEN'S FUNCTIONS 



In this Appendix, we will give expl icit expressions for G{0; Eq) and G(2; Eq) in different energy regions, calculated 
according to the definition Eq. ( 2.35 ). Here we use the same notation as in the text; x = Eq — Eq. When Eq > Ef, 
i.e., X < — 4J(1 + (5^), the energy is above the top of the continuum band, and the Green's functions are 



G{0;Eo) 



^a;2-4J(l + (52)|x| 



(Al) 



G{2;Eo) 



2J(l + 52) l2J{l + S^) 



- 1 



y/x^ -AJ{l + S^)\x\ 



(A2) 



When Ei < Eq < Ef, i.e., — 4J(1 + 5'^) < x < 0, the energy is within the continuum band. 



GiO;Eo) 



^-x2-4J(l + (52)a 



G(2;i?o) 



+ i 



2J(H-52) L2J(l + ,52) J V-a;2-4J(l + ,52) 



(A3) 
(A4) 
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When Eq < Ei, i.e., a; > 0, the energy is in the gap, 

G(0;So) = - 



(A5) 



G{2;Eo) 



+ 



2J(1 + J2) L2J(l + <52) J + 



(A6) 



APPENDIX B: DIAGONALIZATION OF HAMILTONIAN (4.4) 



Hamiltonian (4^) can be written as 
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FIG. 1. Wave functions of the Ag and By, bound exciton with U = IQt, V = 0.5t, and 5 = 0.2. The dashed and solid lines 
correspond to the Ag and By, states, respectively. 
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FIG. 2. Energy of excitons as a function of center-of-mass momentum K with U = 5t, V = 2t, and 5 = 0.2. 



FIG. 3. Optical absorption spectrum with U = 5t and V = t. The Bu exciton occupies 0.520 of the total oscillator strength. 
The arrow indicates the energy of the corresponding Ag exciton. 

FIG. 4. Optical absorption spectrum with U = 5t and V = 2t. The arrow indicates the position of the Ag exciton, and the 
Bu exciton takes 0.949 of the total transition strength. 



FIG. 5. Electron-hole correlation function as a function of the impurity strength Vb for impurity potential i/J™^ with U = lOt 
and V = t. The impurity is situated at site 5 in a iV = 10 site chain. The solid, long dashed, and short dashed lines correspond 
to 7i(5, 6), n{7, 8), and 7^(3, 6), respectively. 



FIG. 6. The hole density at the impurity site vs impurity strength Vb for the impurity potential W^^. The parameters here 
are the same as in Fig. 5. 

FIG. 7. Distribution of charge density in real space for the first type of impurity, f/}™^. The solid circles describe the hole 
density, while the open ones are for the electron density. Here U = lOt, V = t, and Vb = 0.5t. 

FIG. 8. Electron-hole correlation function, as a function of the impurity strength Vo for the impurity potential Z/™^ with 
U = lot and V = t. The impurity is located at site 5 in a A?" = 10 chain. The solid, dashed, and dot-dashed lines illustrate 
7^(5,6), 7^(3,4), and 7^(3,6), respectively. 



FIG. 9. Electron-hole correlation function as a function of the impurity strength Vb for the third kind of impurity J?™^ in 
a, N — 10 chain. Parameters are the same as for Fig. 8, U = lOt and V = t. The solid, dashed, and dot-dashed lines are for 
7l{5, 6), TZ{3, 4), and TZ{3, 6), respectively. 

FIG. 10. Distribution \Zk\'^ for the impurity potential H'^^ with U = lOt and V = t in a N = 12 chain. There are two 
exciton branches {Ag and Bu)- The solid symbols indicate the amplitude of excitons in the Ag branch, and the open ones give 
the amplitude of excitons in the B^ branch. The circle corresponds to Vb = O.lt, the situation before the crossover, and the 
box corresponds to Vb = t, which is larger than the crossover value. 

FIG. 11. Distribution \Zk\^ for the second kind of impurity potential H™'' with U = lOt and V = t in a, N — 12 chain. The 
solid and open symbols have the same meanings as in Fig. 10. The circle describes the case of Vb = O.lt and the box describes 
Vb = t. The line is for the case Vb = 2t. The vanishing amplitude in every exciton state shows the breakdown of the exciton. 

FIG. 12. Distribution IZxf for the third kind of impurity H'^^ with U = lOt and V = t in a N = 12 chain. The circle and 
box correspond to Vb = O.lt and Vb = t, respectively. 

FIG. 13. Electron-hole correlation function as a function of the impurity strength Vb for a charged impurity 7?™^ in a A'^ = 10 
chain. Parameters are the same as for Fig. 8, U = lOt and V — t. The solid and dashed lines are for TZ{5,6) and 7?,(3, 6), 
respectively. 

FIG. 14. Distribution of charge density in real space for a charged impurity with extended potential H'^^. The solid symbols 
describe the hole density and the open ones are for the electron density. The circles and triangles correspond to Vb = 0.5t and 
Vb = 5t, respectively. Here U = lOt, V = t. 
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FIG. 15. Distribution \Zk\'^ for a charged impurity with extended potential if™'' for U = lOt and V = t in a N = 12 chain. 
The soUd and open symbols have the same meaning as in Fig. 10. The circle describes the Vo = O.lt case, the box describes 
Vo = t, and the line is for Vo = 5t. The vanishing amplitude in every exciton state indicates the dissociation of the exciton. 

FIG. 16. Intrachain and interchain exciton bands with U = lOt, V — t, t± — 0.5t and V± — t. The two solid lines describe 
the dispersion of the interchain exciton, and the four dashed lines are for intrachain excitons. 

FIG. 17. Intrachain and interchain exciton bands with U = Wt, V = t, t± = 0.2t and V± = 1.5t. Solid lines are the energy 
spectra for interchain excitons. Here the interchain exciton is the lowest excited state. 

FIG. 18. Wavefunctions of the lower static interchain exciton with U = Wt V = t, t± = 0.2t, and V± = t. The solid line 
gives the relative amplitude that the electron and hole are in the same chain, and the dashed line gives the amplitude that they 
are in different chains. 

FIG. 19. Wavefunctions of the higher static interchain exciton with U = Wt, V = t, t± = 0.2t, and V± = t. Here the 
intrachain amplitude is zero. 

FIG. 20. Intra/interchain electron-hole correlation vs r = V±/V in a two-chain system. Each chain has N = 12 sites and 
the parameters are U = Wt, V = t, and t± = 0.2t. The intrachain correlation function illustrated here by the dashed line is 
7l{5, 6), and the interchain correlation function given by the solid line is TZ{5, 17), where site 12 + i indicates site i in the second 
chain. 

FIG. 21. Schematic diagram of the charge transfer from a polymer chain to a dopant molecule. The dopant has an acceptor 
level which interacts with the polymer by nearest-neighbor hopping v. 

FIG. 22. Probability of charge transfer of the Bu state as a function of acceptor level Ae. The parameters in the polymer 
are U = Wt, V = t, and 6 = 0.2, corresponding to the energy of the B„ exciton, Eo = 9.537t (indicated by the arrow). The 
dopant concentration c = 0.2. In the solid line, v = 0.212t, and in the dashed line, v = 0.566t. 

FIG. 23. Probability of charge transfer of the Ag state as a function of acceptor level Ae. The parameters are the same as 
those in Fig. 22, and the corresponding Ag state has energy Eq = 9.140t (indicated by the arrow). The solid and dashed lines 
correspond to v = 0.212t and v = 0.566i, respectively. 



24 



1.2 



0.6 - 



0.0 • ® • o • % % • 



o o 

-0.6 - 



-I 2 ' ^ ^ ^ ^ ^ ^ ^ ^ ' 



0.6 



0.3 



N 



0.0 



8 • 



o 



□ 



-0.3 



-0.6 



-1.0 



-0.5 



0.0 
K/Tl 



0.5 



1.0 



0.6 



0.3 



N 



0.0 



-0.3 



-0.6 




-1.0 



-0.5 



0.0 
K/Tl 



0.5 



1.0 



0.6 



0.3 



N 



0.0 



-0.3 



-0.6 



-1.0 



m • m 



o 



-0.5 



0.0 
K/Tl 



0.5 



1.0 



1.0 



0.5 



0.0 • i * • 



-0.5 



-1.0 ' ^ 

12 3 4 



I 



8 



A 



A 



7 8 9 10 



0.6 



0.3 



OA 



N 



0.0 



-0.3 



-0.6 



-e- 



-1.0 



-0.5 



0.0 
K/TT 



0.5 



1.0 



10.0 




8.5 



8.0 

0.0 



0.5 
K/71 



1.0 



10.0 



9.5 



9.0 





8.0 

0.0 



0.5 
K/71 



1. 



-OA 



-10 



-5 




I 



5 



10 



1.0 



0.5 



0.0 



0.0 



0.5 



1.0 



1.5 



polymer 



dopant 



